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1. Introduction 

Topological properties of QCD are among those most widely studied on the lattice. Various 
methods have been used to this end, involving renormalization, cooling, fermionic zero modes, 
geometric definitions, etc. In recent years, the advent of fermionic actions, such as the overlap, 
which do not violate chirality, has brought a new thrust to the subject. 

In this work we compute the renormalization constants which are necessary in order to ex- 
tract topological properties, in the "field theoretic" approach, from Monte Carlo simulations using 
Wilson or Symanzik improved gluons, and clover or overlap fermions. We compute the multiplica- 
tive renormalization Zg of the topological charge density, to 1 loop in perturbation theory, and the 
power divergent additive renormalization M(g 2 ) of the topological susceptibility, to 3 loops. 

The main motivations for doing this work are: a) To enable comparison between different 
approaches used in studying topology, so that a coherent picture of topology in QCD may emerge, 
b) To enable studies, in numerical simulations, of quantities involving the density of topological 
charge, q[x), rather than only the integrated charge; this is necessary, e.g., for studying the spin 
content of nucleons. c) As a feasibility study in lattice perturbation theory: Indeed, this is the first 
3-loop calculation to appear in the literature, involving overlap fermions. 

A more extensive write-up of this work can be found in our Ref. [JXJ] . 



2. Computation of Zg 

Our first task is to compute the multiplicative renormalization Zq [Q] of the topological charge 
density q L (x) to one loop, using the background field method. We use the standard definition of q L , 
given by a product of twisted plaquettes 

1 ±4 

^l( x ) = -^=2 E e M vpffTr[n M v(x)n p(T (x)] (2.1) 
z 71 }ivpo=±\ 

(£-Ai.v,p,ff = — e fi,v,p,<r ; in standard notation: Yl llv (x) = U fl (x)U v (x + n)Uj l (x + v)Ul(x).) 

The classical limit of the operator shown in Eq.(2T) must be corrected by including a renor- 
malization function Zq, which can be expressed perturbatively as 

Z Q = l+Z l -g 2 + --- , Z 1 =Z n -N c +Z l2 /N c (2.2) 

We perform a calculation of Z\ ; this involves only the gluon part of the action. 
In the background field method, link variables are decomposed as 

U^x) =V^x)U c ^(x), V tl (x)=e* a M, U cll {x)=e iaB ^ (2.3) 

in terms of links for a quantum field and a classical background field, respectively. 
The diagrams involved in the one-loop calculation of Zq are shown in Figure 1 . 
The standard Symanzik improved gauge field action, with 4- and 6-link Wilson loops, is 

2 r 

Sg=— c J- R&Tr ( l -^plaquette) + C\ £ Re Tr ( 1 ~ ^rectangle ) 
& plaquette rectangle 

+ C 2 Re Tr ( 1 - ^chair ) + C 3 Re Tr ( 1 - i/parallelogram ) (2.4) 

chair parallelogram 



230/2 



Topological charge renormalization 



Apostolos Skouroupathis 






(a) (b) (c) (d) 

Figure 1: Diagrams contributing to Z\ . The bullet stands for topological charge vertices. 



The Symanzik coefficients c,- must satisfy: cq + 8c i + I6C2 + 803 = 1. 

Our calculations are performed without assumptions on the values of the external momenta 
Pi, P2 '■ This is safest for q L , otherwise one may easily end up with indeterminate expressions; it 
entails handling 3-point form factors ^ in D = 4 — 2e dimensions, such as: 

(ka) 2e f d D k sink^ sink v 

C, v (a, PuP2 ) = —J (27r) ^ 2( ^ i)2( , + ^^ 2)2 ( 2 - 5 ) 

Diagrams (c) and (d) of Figure 1, taken separately, exhibit poles in e ( (d) oc — \/e — l n K 2 a 2 ). 
These cancel, however, upon summation, as is expected by the fact that Q does not renormalize 
in the continuum. The calculation of Zq is particularly involved in the present case, involving 
propagators and vertices from the improved gluonic action. In particular, the calculation of diagram 
(d) involves a summation of > 1000000 different algebraic expressions at intermediate stages. 

Our results for Zq are listed in Table 1. In all calculations that involve the parameters c,, we 
choose a standard set of values, as in Ref. [Q]. The choice of the sets of parameters correspond to the 
most popular actions: The first set corresponds to the plaquette action, the second set corresponds to 
the tree-level Symanzik improved action [Q] and the next 6 sets correspond to the tadpole improved 
Luscher-Weisz (TILW) action [§] for 6 values of beta: j8 = 8.60, 8.45, 8.30, 8.20, 8.10, 8.00 . The 
last two sets correspond to the Iwasaki [f7|] and DBW2 ^ actions, respectively. 

In the case of the plaquette action, our result agrees with the known result of Ref. [g]. 

It is worth noting that the value of Z\ (and of e$, see below) for the DBW2 action is the 



smallest one, leading to a renormalization factor Zq closer to 1 (and M(g ), Eq. ( |3.2| ), closer to 0). 



This would single out the DBW2 action as a better candidate for studies of topology. 

3. Computation of M(g 2 ) 

The second task we attend to is the calculation of the additive renormalization of the topolog- 
ical charge susceptibility, which is defined as 

X L = (0)> (3-1) 

X 

X L develops an unphysical background term which becomes dominant in the continuum limit 

X L (g 2 ) = a%(g 2 ) 2 x+M(g 2 ), M (g 2 )=e 3 -g 6 + e 4 -g g + --- (3.2) 
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Figure 2: 2-loop diagram contributing to <?3. Bullets stand for topological charge vertices. 





(a) 



(b) 



Figure 3: Diagrams contributing to e 4 . Wavy (straight) lines correspond to gluons (overlap/clover fermions). 



M(g 2 ) is the power divergent additive renormalization of % L . We compute the 2-loop coefficient 
^3 . This quantity is evaluated for several sets of values of the Symanzik improvement coefficients. 
Figure 2 shows the diagram contributing to . 



The 3 -loop term £4 of the expansion of M(g ) equals e\ 



J 



e 4 + e 4 > where e{ stands for the 



fermionic contribution to e\ (csw is the coefficient in the clover action) 

e{ = N f {N? -l)N c - (04,0 + C4,l csw + ^4,2c| w ) 
and e\ is the purely gluonic contribution, expressed as in Ref. [^] 



(3.3) 



16 



(AT 2 -1)(1.735a£ -10.82 + 73.83//V 2 ) x 1(T 7 



(3.4) 



In fact, what we are interested in, is the calculation of the parameters e^o, , e^%. This task 
is performed using both overlap and clover fermions (clearly, overlap fermions involve only the 
parameter £4.0). Figure 3 shows the 3-loop diagrams contributing to the evaluation of e£. 

The propagator and vertices of overlap fermions can be obtained from the following expression 
for the overlap action, written in terms of the massless Neuberger-Dirac operator [ 10] 



Soverlap = a 4 £ (n)D N (n,m)m {m) 

n,m 



D 



•n 



Mo 
a 



1 + 



X 



Vx^x 



(3.5) 



Mq is a real parameter corresponding to a negative mass term. Mq must lie in the range < Mq < 2r, 
r being the Wilson perameter (in our case r = 1). X is the Wilson-Dirac operator with mass — Mq . 



The clover (SW) fermionic action [|11|], contains an extra term, parameterized by a coefficient, 
csw ; this coefficient is treated here as a free parameter. 

In performing this calculation, a large effort was devoted to the creation of an efficient 3-loop 
"integrator", that is, a metacode for converting lengthy 3-loop integrands into efficient code for 
numerical integration. The output code of the integrator precalculates a number of time-consuming 
common ingredients (Symanzik propagator, overlap expressions, etc.), exploits symmetries of the 
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integration region, integrates in parallel over non-overlapping loops, organizes the integrand as an 
inverse tree for optimized evaluation of innermost loops, etc. 

Table 1 contains our results for (cf. Eq.([3~2])) for different gluonic actions. These results, 
for the case of the plaquette action, agree with older known results (see, e.g., [|l2|]). Figure 4 shows 
the coefficients e^o, 64,1, £4,2 of the clover result for different values of the bare fermion mass m. 
Figure 5 exhibits the dependence of e 4 , using the overlap action, on the parameter Mq. 

A complete tabular version of our results on e\ , for both the clover and overlap cases, can be 
found in our longer write-up, Ref. [jjj]. 





c\ 


C3 


Z n 


Zn 


e-$ x 10~ 7 


1.0 


0.0 


0.0 


-0.33059398205(2) 


0.2500000000(1) 


6.89791329(1) 


1.6666666 


-0.083333 


0.0 


-0.2512236240(1) 


0.183131339233(1) 


3.1814562840(7) 


2.3168064 


-0.151791 


-0.0128098 


-0.20828371039(3) 


0.147519438874(3) 


1.8452250005(2) 


2.3460240 


-0.154846 


-0.0134070 


-0.20674100461(1) 


0.146259768983(1) 


1.8054229585(4) 


2.3869776 


-0.159128 


-0.0142442 


-0.20462181183(1) 


0.144531861677(4) 


1.7516351593(8) 


2.4127840 


-0.161827 


-0.0147710 


-0.20331145580(1) 


0.143464931830(1) 


1.7188880608(5) 


2.4465400 


-0.165353 


-0.0154645 


-0.20162651307(1) 


0.142094444611(2) 


1.6773505020(9) 


2.4891712 


-0.169805 


-0.0163414 


-0.19954339172(1) 


0.140402610424(1) 


1.626880218(1) 


3.648 


-0.331 


0.0 


-0.15392854668(1) 


0.105132852383(2) 


0.752432061(7) 


12.2688 


-1.4086 


0.0 


-0.0617777059(4) 


0.038277296152(6) 


0.04881939(4) 



Table 1: The values of Zn and Z12 (Eq.(2.2), Figure 1), and of ej, (Eq.(3.2), Figure 2), with Symanzik 
improved gluons, for various values of the coefficients cq, c\ , C3 . (C2 = 0) 
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Figure 4: Variation of the terms contributing to e 4 as a function of m 
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Figure 5: Value of e^ as a function Mq 
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